Condensates induced by interband coupling in a double-well lattice 
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We predict novel inter-band physics for bosons in double- well optical lattices(OL). An intrinsic 
coupling between the s and p x band due to interaction gives rise to larger Mott regions on the phase 
diagram at even fillings than the ones at odd fillings. On the other hand, the ground state can form 
various types of condensates, including a mixture of single-particle condensates of both bands, a 
mixture of a single-particle condensate of one band and a pair-condensate of the other band, and 
a pair-condensate composed of one particle from one band and one hole from the other band. The 
predicted phenomena should be observable in current experiments on double- well OL. 



I. INTRODUCTION 

The last decade has seen many exciting developments 
in the studies of cold atoms in OlP. A number of many- 
body phases of bosons and fermions have been observed 
in experiments^^, showing that cold atoms in OL pro- 
vide an ideal platform for studying a number of theoret- 
ical models in condensed matter physics. For example, 
bosons in OL have been used to study the single-band 
Bose-Hubbard model extensively in the past few yeariP^. 

Despite the fact that single-band models are usu- 
ally good descriptions for atoms in OL, recent stud- 
ies have found important influence from the higher 
bands 89 . The leading-order effect of such inter-band 
couping is to renormalize the parameters of the single- 
band modeP^HMl. p or current experimental parameters, 
the phase boundary is found to change slightly while the 
structure of the phase diagram remains the same^l. It 
is therefore commonly assumed that inter-band coupling 
only leads to small quantitative changes in the quantum 
phase diagram of bosons in OL. 

A question then naturally arises, whether OL can be 
used to create many-body systems, which go beyond the 
description of the renormalized single-band models, i.e., 
whether situations exist where interaction induced inter- 
band coupling could give rise to qualitative new quantum 
phases not conceivable in the single-band picture? It is 
easy to see that the chemical potential \i needs to be much 
smaller than the band gap A^ between the s and p band 
for invalidating the standard single s band description. 
In this Letter, we propose using a system in which A^ 
can be tuned to be comparable with or even smaller than 
\i. A single-band description for this system is no longer 
valid, and novel quantum many-body phenomena due to 
inter-band coupling emerge. 

The system we consider is bosons in a double-well lat- 
tice. In the literature, there has been studies on bos ons in 
an excited band in double-well OL or standard Oil 14 * 15 *. 
Here we focus on a new aspect regarding the inter-band 
physics between the lowest two bands in a double-well 



lattice. From the experimental side, double- well lattices 
have recently been realized by several groups with dif- 
ferent lattice geometries'^. We focus on the simplest 
case, where a second laser is added to a simple cubic 
optical lattice along the x direction 18 . The physics we 
discuss here can be easily generalized to other cases. By 
choosing the wave length of the second laser as half of 
the first one which creates the cubic lattice, a superlat- 
tice forms with a double-well potential. For appropriate 
choices of the relative positions of the two lattices, the 
potential is symmetric with respect to the center of each 
lattice site, as shown in Fig.([]Ji). 



II. HAMILTONIAN 

The potential of a double-well lattice in our case 
can be written as V(R) = V X (R X ) + V y (R y ) + V Z (R Z ), 
where V X (R X ) = Vl sin 2 ( ZL ^ 3L ) — Vs sin 2 ( ) , V y (R y ) = 
^sin 2 (^), and V Z (R Z ) = ^sin 2 (^f-), V L (ft) and V s 
are the amplitude of the long and short lattice along the 
x(i = y, z) direction, d is the lattice spacing. For this 
periodic potential, the single particle energy e cr (k) and 
the Wannier wave function i/j a (R) can be exactly calcu- 
lated, where a is the band index for the three dimensional 
lattice. e a (k) can be written as the sum of the corre- 
sponding single particle energy e^ n (A^), where i = x,y,z 
and n is the band index, for the one dimensional lattice 
Vi(Ri). Typical band structures for a standard optical 
lattice and a double- well lattice are shown in Fig.Q. It 
is clear that in a double- well lattice, e x ^(k x ) and e Xj i(k x ) 
become very close to each other and are well separated 
from higher ones. Thus we focus on the lowest two bands 
for the three dimensional lattice, the energy of which 
can be written as e s (k) = e x $(k x ) + e y $(k y ) + e z $(k z ), 
e Px (k) = e Xi i(k x ) + e yt o(k y ) + e z ,o(k z ). 

The Wannier wave function for each band can 
be written as ip s (R) = ^x$(Rx)^y$(R y )^z$(Rz), 

tl>Px(R) = ^x,i(Rx)^y$(R y )^z$(Rz), where ^ n (Ri) 
is the Wannier wave function for the one dimensional 
lattice Vi(Ri) as shown in Figffib). Because the 
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FIG. 1: (a) Schematic of a double- well lattice along the x 
direction, as well as the lowest a few bands marked down in 
different colors. Contours with an unit cell marked in a red 
box on the x — y plane are also shown. Dark region represents 
the potential minimum. (b) Wannier functions for the one di- 
mensional lattice Vx(Rx), where Vs — VLEr, Vl = 4Er, V y = 
V z = 13Er along the x direction. 



lattice potential V(R) respects the inversion symme- 
try, i.e., V(R) remains unchanged under the transfor- 
mation Ri — > — Ri, ip a (R) has well defined parity. 
For example, ip s (R x , R y , R z ) = ip s (-R x , R y , R z ) and 
ip p (R x , Ry> Rz) — —^p(—Rx,Ry,Rz)' Both of them re- 
main unchanged for R y — >> —R y or R z — >• —R z . This 
is similar to the parity of the usual s and p x orbitals in 
a single harmonic oscillator. This is the reason that we 
denote the lowest two bands as the s and p x bands for 
simplicity. 

The many-body Hamiltonian can be written as H = 



m ' k.sm ' °pm 



H C1 where 



amr cm 

is the usual Bose-Hubbard band model for the s and p x 
band, and 



E 



(wh; 



C.C^ 



(2) 



is the inter-band paring coupling term. In above two 
expressions, the band index p x has been replaced by p 
for simplicity, r — dx,dy,dz represents the unit vec- 
tor along the x, y, z directions, m creates a particle 
in the basis of Wannier wave functions at a lattice site 
m = d(m x , rriy , m z ) that contains two wells, \i the chem- 
ical potential, fi s = /i, \i v = fi — A g: A g is the energy 
difference between the middle of the s and the p x band, 
which will be referred to as the band gap for simplicity. 
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FIG. 2: (a) The energy spectrum for a standard lattice, where 
V L = 13E R , V s = 0, = 6E R . (b) The energy spectrum for 
a double-well lattice, where Vl = 4Er,Vs = 11Er,A 9 = 
1.5Er in ID along the x direction. For comparison, the dif- 
ferences between the maximum and the minimum of these 
two lattices have been chosen to be the same. Red and blue 
horizontal lines in the upper panel represent the energy level 
of each Wanner orbitals. 



t a f = J2k ^a(k)e lk ' r /N s is the tunneling amplitude be- 
tween the nearest neighbor sites, where N s is the total 
number of lattice sites. From exact numerical results, we 
found that tunneling amplitude beyond nearest neighbor 
sites is at least two orders of magnitude smaller, and thus 
can be ignored(See Appendix). 

The interaction part in the Hamiltonian was ob- 
tained by expanding the field operator ^(R) in 
the general expression for the interaction of bosons 
U = / d 3 R^(RY^(R)^(R)^(R) using 4f(R) = 

^2 am ^a(R — m)^crm- A straightforward calcula- 
tion shows that there are three terms for on- 
site interaction. ^Z am nam(n am — 1) is the 
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FIG. 3: Fock states and the matrix elements between them 
due to on-site interaction for two particles per site. 
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FIG. 5: Mean- field phase diagram, where Vs = HEr, Vl = 
4Er, V y = V z = 25Er. Phase boundary between a conden- 
sate phase and the Mott insulator is shown in different colors 
for different filling factors in the Mott insulator. 



FIG. 4: (a) A g (purple dash-dotted), tunneling along the x di- 
rection of the s (red solid) and p x band (blue dashed) as func- 
tions of Vs, where Vs + Vl = 15Er. (b) interaction C7 s (red 
solid), U p (b\ue dashed), U sp (brown dotted) and W (purple 
dash-dotted), where a s /d = 0.05 and V y = V z = 25Er. 



intra-band repulsion for each band, where U a = 
f d 3 R\ilj a (R)\ 4 . The second one is the inter- 
band repulsion, U sp = U sp ^ m h sm h pm , where U sp = 

^^fd 3 ii|^(jR)| 2 |^p(jB)| 2 . The last one as shown 
in Eq.|2| describes the interacting process that scatters 
two atoms from the s band to the p x band, where W = 
^> / d s R^f(R)^l(R). Because of different parities 
of the s andp x band, ip s (R x ,R y , Rz) = i(j s (—R x , R yi R z ), 
ijjp(R x , R y , R z ) — —ipp(—Rx,Ry,Rz), the Hamiltonian 
does not contain terms like n a b^b a ^ where a ^ a' . We 
have dropped off other terms involving inter-site inter- 
action that are much smaller. The Fock states for two 
particles per site as well as the matrix elements between 
them due to on-site interaction are shown in Fig. ([3]). 

The Hamiltonian in Eq.Q has earlier been considered 
for standard OL with a fixed band gap at a constant 
lattice depth in a context different from our interestP^HHl 
Dependences of the parameters of the Hamiltonian in 
Eq.Q on Vl/E r and Vs/E R for a typical double- well 
lattice are shown in Fig.( ~~ 
the recoil energy. It is c 
decreases quickly. 



4j), where E R = h 2 /8Md 2 is 
ear that as Vs increases, A^ 



III. MEAN-FIELD PHASE DIAGRAM 

To solve this problem, the Gutzwiller mean field ap- 
proach is used. The trial wave function is written as 
hA) = rimZ)a,/3 c a/3l a ^)m, where |a,/3) m is the on-site 



Fock state, a and f3 are the particle numbers in the s 
and p x orbital respectively. The order parameter for each 
band (b am ) = (j) a is determined self-consistently, and the 
phase boundary is obtained from \(j) s \ = \<fr p \ =0. 

The phase diagram for the parameters Vs = H£^r, 
V L = AE R and V y = V z = 2bE R is shown in Fig.([5|). 
Compared with the usual results for the single-band 
Bose-Hubbard model, it is clear that the structure of the 
phase diagram changes dramatically. For the single-band 
model, the area of the Mott region on the phase dia- 
gram monotonically shrinks with increasing filling factor. 
For our two-band model, there is an "even-odd" effect, 
namely, the Mott region at an even filling 2uq becomes 
larger than the one at an odd filling 2no — 1. 

This "even-odd" effect can be first understood intu- 
itively using an extreme case, Vl = and Vs ^ 0. For 
this particular configuration, the Mott regions with odd 
fillings must vanish, since the lattice spacing has changed 
to d/2 and odd fillings in the original basis now corre- 
spond to non-integer fillings. However, for any general 
case Vl ^ 0, the lattice spacing is still d. The above sim- 
ple picture does not apply. To fully understand the "odd- 
even" effect, one needs to consider on-site correlations 
induced by the inter-band coupling and the competition 
between the onsite energy and the tunneling energy. 

The excitation gap of a Mott insulator can be writ- 
ten as A n = E n +i + E n _i — 2E ni where E n is the on- 
site energy of the Fock state with n particles per site. 
In a single s-band model, E® = U s n(n — l)/2 — fin. 
Consequently, A° = U s ~ a s does not dependent on 
n, where the superscript represents the results for the 
single band model. On the other hand, the tunneling 
energy — £ s (6| m 6 sm +i) ~ —t s ^/Jn^rY)n increases with 
increasing n due to bosonic enhancement. As a result, 
the critical value of a S: i.e., at the tip of the Mott re- 
gion, has to increase so as to overcome the increasing 
kinetic energy with increasing n and form an insulator. 
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This is the reason that Mott region monotonically shrinks 
on the phase diagram of the single band model with in- 
creasing n or \i. In our two-band model for double- well 
lattices, H c lowers the onsite energy E n by coupling the 
s and p x band, E n = E® + 5E n , and modulate the na- 
ture of the excitation gap A n = A° + 5A n . For ex- 
ample, SE 2 = -4W 2 /(2A g + U p - U s ) + 0(W 4 ) and 
SE 3 = -12W 2 /(2A g ^ + U p + 2Z7 sp - 3/7 s ) + 0(W A ) be- 
cause of the formation of the pair C2o|2,0) + co2|0,2) 
and C3o|3,0) +ci2|l,2) respectively, where 020(^30) and 
£02(^12 ) have a phase difference of ir. SE 2 (SEs) can be 
viewed as an effective binding energy induced by the par- 
ing coupling H c with respect to |2, 0)(|3, 0)). On the 
other hand, 5E\ = as the state |1,0) or |0, 1) does not 
couple to any other Fock states. One can easily verify 
that SAi < SA 2 if U sp > U a , where SA 1 = SE 2 < 0, 
SA 2 = 5E 3 - 25 E 2 . Thus A 1 = U s + SA 1 is suppressed 
and meanwhile A 2 = U s -\- 5A 2 is either enhanced or re- 
duced less than Ai, i.e., Ai < A2. As seen from Fig. ([4] 
(b)), U sp > U s is satisfied if Vs > 2.5Er. Moreover, 
U sp — U s increases with increasing Vs, so does A2 — Ai. 
Though populating the p x band may increase the kinetic 
energy for n = 2 compared with n = 1, when Ai < A2, 
this kinetic energy increase, as well as the bosonic en- 
hancement, is no longer important, and the critical value 
of a s for n = 2 Mott insulator becomes even smaller than 
then one for n = 1. That is what we see on Fig. (JSJ). Sim- 
ilar analysis applies for n = 3, 4, and so on. 

This "even-odd" effect can be directly observed from 
in-situ images of the density of trapped bosons in a 
double-well lattice. With increasing interaction, the 
Mott plateau with two atoms per site would first emerge 
when Vs is sufficiently high. The above discussion can be 
generalized to arbitrary orbital numbers n° in each lat- 
tice site. When there are strong pairing couplings similar 
to H c between different orbit als, the Mott region on the 
phase diagram for (n mod n°) = will be larger than the 
one for (n mod n°) = 1, 2, .., n° — 1. Similar phenomenon 
also happens for a model containing two or three p bands 
in standard OlPEl. 



IV. CONDENSATE PHASES 




FIG. 6: (a) n s + n p (brown dash-dotted), n s (red solid) and 
n p (blue dashed) versus \i at a s /d = 0.081 for Fig.|5}. 
(b) (b s )(ied solid), (& s fr s )(pink dash-dotted), (b p b p ) (brown 
dashed) and (b p ) — for same parameters as (a). 



A. {CI) condensate 

For a (CI) condensate to emerge, the chemical poten- 
tial fi needs to be larger than a critical value /i c . We 
write down the energy density for the excited band in the 
mean field approach as E p = —(/J* — A g + 2 t p ?) \ (j) p \ 2 + 
U sp n s \(l)p\ 2 — 2W\(j) s \ 2 \4>p\ 2 + U p \(j)p\ A . Based on Landau's 
criteria, fi c = A g — 2^^t p f + U sp n s — 2W\(j) s \ 2 and 
\(j) p \ ^ for /i > p c . In the parameter regime consid- 
ered in Fig.Q, the very large U sp has largely increased 
the value of of \i c . Thus, the ground state forms a (C2) 
condensate with (b sm ) 7^ 0, (b pm ) = and (bp m bp m ) =^ 0. 
At lower lattice depths, i.e., Vl = 4Er and Vs = 1 
U sp is small enough and (CI) condensates can emerge, 
i.e., for about four atoms per site 26 . 



Now we turn to the condensate phases. Because of 
more than one orbital in each single site, there are mul- 
tiple choices for the bosons to form condensate phases. 
(1) When |(6 sm )| ^ and |(6 pm )| + 0, we call it a (CI) 
condensate, which is a coherent mixture of single-particle 
condensates of both bands. (2) When one of the single- 
particle condensates vanishes, we call it a (C2) conden- 
sate, i.e., |(6o- m )| = an d |(&<7'm)| 7^ 0> we cau it a (C2) 
condensate. (3). When |(6 sm )| = |(&pm)| = 0, the inter- 
band coupling induces a paired-condensate \ (bl m b pm } \ = 
|(frj m 6 sm )| 7^ under certain conditions. 



B. (C2) condensate 

In double- well lattices, a (C2) condensate has a unique 
property: when (b am ) 7^ 0, (frcr'mfrcr'm) must be finite(see 
Fig. ^). We can view the (C2) condensate as a coher- 
ent mixture of a single-particle condensate of one band 
and a two-particle condensate in the other band. A non- 
zero (bp m bp m ) despite a vanishing (b pm ) results from the 
inter-band coupling H c . Once (b sm ) =^ 0, the ground 
state includes Fock states \a, /3) m and |a±z/, /3) m , where 
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v is an integer. On the other hand, the coupling term 
H c provides a matrix element between \a ± 2z/, /3) m and 
|a, /? =F 2z/) m . The expectation value for 6 pm 6 pm then 
must be finite, since the ground state includes both 
\a,/3) m and |a,/3 ± 2^) m , though it does not include 
Fock states |a,/3 ± (2v + l)) m , i.e., (6 p ) = 0. This is 
very different from a single-band model, where (b) and 
(66) must vanish at the same time. Here, (6 cr / m 6 -/ m ) can 
be induced by a finite (6 crm ) in a different band. For 
standard OL, a (C2) condensate also exists but the am- 
plitude of (b pm b pm } is rather small because of the very 
large band gap. Our double- well lattices significantly re- 
duce the band gap and make ( C2) condensate observable 
in experiments. 

A (C2) condensate can be probed in Time- 
Of-Flight(TOF) experiments with a band mapping 
technique 27 , which maps the crystal momentum distri- 
butions in different bands, i.e., (n s ^) and (n p ,k), sepa- 
rately to the first and second Brillouin Zone after expan- 
sion. Since there is no single-particle condensate in the 
p band, i.e., (6 pm ) = 0, no singular feature is present 
in (%>,k)- On the other hand, a finite (b pm b pm ) induces 
a novel structure in the momentum correlation function 
ip(ki,k 2 ) = (np, kl n p , k2 ) - (h pMl )(h pM ). A straight- 
forward calculation shows, / p (ki,k 2 ) = I\ + I 2 , where 
h = IE mi (n mi )e i(kl - k2) ' mi | 2 - E mi (^ mi ) 2 and I 2 = 

IE mi (^ mi ^ mi )e i(kl+k2) ' mi | 2 " E mi l(^ mi ^, mi )| 2 - 
As a function of ki and k 2 , h has singular peaks at 
ki + k 2 = <?K, where q is an integer, and K = 2 J L (^, y, z) 
is the reciprocal lattice vector. By measuring the ampli- 
tude of I p (Q,— Q), where Q is an arbitary momentum 
in the second Brillouin Zone, the strength of (6£ mi &£ mi ) 
can be probed. 

It is worthwhile to mention that, to observe (C2) con- 
densates discussed here, the temperature of the atoms 
needs to satisfy T < t s ,W. For a typical recoil energy 
Er ~ lOOnif , it can be seen that t s and W are of the or- 
der of tens nK. It is therefore quite promising to observe 
(C2) condensates and the novel structure of the phase 
diagram in current experiments. 



C. (C3) condensate 

The Fock space in each single site can be divided 
to two subspaces, {Si} : \n — 2/,2/) m and {S 2 } ' 
\n — 21 — 1, 21 + l) m , where I = 0, 1, 2, since H c only 
mixes states within each subspace. For n particles per 
site, we can denote the ground state of each subspace 
as |^ n i)m and |^ n 2)m and the ground state energy as 
E n i and E n2 . Both E n \ and E n 2 are dependent on 
a s and the lattice configuration. For the parameters in 
Fig.(§, E n 1 <C E n2 . On the other hand, under certain 
conditions, E n i « E n 2. For example, for Vs = 1£^r, 
Vl = V y = V z = 80Er, which is essentially an one 

dimensional system, E n \ « E n 2 is satisfied for n = 3(5) 




FIG. 7: Schematic of a second order process that induces the 
correlations between the nearest neighbor sites. The bounded 
dimers represent the state C2o|2, 0) +C02IO, 2). The intermedi- 
ate state with one and three particles in the nearest neighbor 
sites has much higher energy. The red vertical arrows repre- 
sent the pseudospin in the spin model. 

around a s /d = 0.18(0.123). Under this situation, by tak- 
ing into account higher order correlations between near- 
est neighboring sites, an exotic (C3) condensate arises. 

We apply a spin-1/2 representation for the system, 
I t)m = |^ni)m and I l) m = |^n2)m. The spin oper- 
ators are defined as S^ = |(| t)mm(t I - I l)mm(l|), 
= I t>mm(l |, S~ = | |>mm(t I, which sat- 
isfy the commutation relations, [S+,S~] = 2S^ n £ mjn , 

^n] = ^m^m.n- [S m , S„] = ^ m ^m.n- A Spin-l/2 

model can be obtained by applying second order pertur- 
bation theory(see Fig.|7|), H spin = Y, m ,r J A S m S m +? + 

c ' c ) + Xlm,r ^ri^mSm+r + C ' C ) + Em,f ^^m^m+r — 

/iEm^ where J^A^ - t s ft p f/(E ex - E nl - E n2 ), 

J Z r ~ tlr/i^ ~ E nl ~ E n2 ) + ~ E nl ~ E n2 ) , 

h « E n2 — E n i and E ex is the energy of an excited state 
with n — 1 and n + 1 particles on nearest neighbor sites. 

Unlike the usual XX Z model obtained from a single- 
band Bose-Hubbard model, there is an additional term 
(^m^m+r + C ' C ) m our ca se, since the pseudospin is rep- 
resented by Fock states with the same particle number 
and thus the total spin S z = ^ S^ does not need to 
be conserved. Because of this additional term, the spin 
model is no longer invariant under U(l) gauge transfor- 
mation, i.e., S+ — ^ S+e* 7 where 7 is an arbitrary phase. 
Instead, only at discrete values of 7, i.e., 7 = ±7r, the 
spin model is invariant. This indicates a Z 2 symmetry of 
this spin model. 

A mean field approach can be used to solve the above 
spin model. For example, when JV, J5 < 0, written down 
{SD = l cos0 ' (^m) = isinfle 4 * (S~) = Isinfle"^, 
the mean field energy of the spin model can be expressed 
as (H spin ) M = sin 2 0(J t + A t cos 2cp)/2 + ( J t z cos 2 (9)/4 - 
(hcos0)/2, where J t = Jr, A t = J2r^r and Jf = 
J*. It can be easily verified that when \J t \ + |A^| > 
\Jf\/2 is satisfied, there is a critical value h c = \Jt\ + 
|A t | — \ J{;\/2. If \h\ < h C: the value of to minimize 
(H spin )M is no longer zero or 7r, and there is a projection 
of the spin on the x — y plane with a preference on the 
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direction, i.e., tp = ±tt/2 for A t > or cp = 0, tt for 
At < 0. Similar results can also be obtained for other 
signs of J>, J5. 

We now go back to the boson representation 

\G) = n m ((^4) i l*ni)m + (sinf)3e^|* n2 ) m ). As 

|^ni)m = ^2iCu\n-2l,2l) m and |^ n2 ) m = Z)jC 2 z|n- 
21 — 1,2/ + l) m , it can be seen that (5+) ^ implies 

(^m^p.m) = (>,m^s,m) 7^ 0- We refer to this state as a 
(C3) condensate that emerges when \h\ < h c . It is a pair- 
condensate composed of one particle from one band and 
one hole from the other band. (C3) condensates have also 
been found in a model describing two component bosons 
in a lattice, however, due to different underlying physics 
that belong to different universality class, i.e., breaking 
U(l) other than Z2 symmetry as in our case^. We would 
like to point out that it is challenging to observe a (C3) 
condensate in current cold atom experiments as it re- 
quires reaching a very low temperature, T ~ t a t a f/U a , 
which is of the order of 0.1 — InK. 

Finally, the above spin model in ID can be mapped 
to a fermion model via Jordan- Wigner transformation, 

H F = E(mn) J ^(c]nC n +C.c) + E (mn) A mn (c^gt + C.c) + 
4 ^E(mn)4cm4cn ~ (2ft + 4JJ~) ^ m C^m, where 

A mn = — A mn = Adx- The fermion operators satisfy 

Cm = (n n <m S n) S L and 4i = {U n <m S n) 5 ~ • Except 

for an unimportant nearest neighbor interaction term, 
this fermion model is identical to a model describing 
one dimensional p-wave superconductor 29 . The degen- 
erate ground states in the fermion model correspond to 
(p = for At > or 0/ir for A t < in the boson 

model. 



V. CONCLUSIONS 

We have shown that the band gap between the s and 
one of the p bands can be easily tuned in a double-well 
lattice. This tunable band gap gives rise to novel 
properties for bosons. Interaction induced inter-band 
coupling becomes crucial in this system and modifies 
drastically the nature of the phase diagram, as well as 
the possible condensate phases. For the phase diagram, 
we found an "even-odd" effect, namely Mott region with 
a filling factor 2no on the phase diagram becomes larger 
than the one with an odd filling 2no — 1. We also found 
three different types of condensate phases. Particularly, 
when strong interaction suppresses the fluctuations of 
the total density, inter-band coupling can nevertheless 
induce an paired-condensate that is composed by one 
particle from one band and one hole from another band. 
Our results indicate the strong possibility of realizing 
novel condensate phases in double- well OL. Our studies 
can be generalized to the cases that more than one 



higher band is tuned to be very closed to the lowest 
band, where very rich physics is expected. 
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VII. APPENDIX 

The kinetic energy in the momentum space can be 
written as JC = Y] re zc^ ^d r, where a) -(a A cre- 

t-^ok ak a k ak 1 ak' 

ates (annihilates) a particle in the basis of Bloch wave 
functions </> ai (R). Define ftt m = £g ^g^/V^VT, 
where N s is the total number of lattice sites, which cre- 
ates a particle in the basis of Wannier functions VVm(^) 
at site m, JC = Eamri V/^mW + c.c), where / is 
an positive integer, r = dx, dy, dz represents the unit vec- 
tor along the x, z directions, and tunneling constants 

ta,r,i = Sit e a k e%lk ' r /N s - Because the single particle en- 
ergy is already diagonal with band indices, there is no 
crossing term involving tunneling between orbitals with 
different index a. It can also be shown from a straight- 
forward calculation that t a f,i corresponds to the overlap 
integral of the Wannier wave functions ?/v (R — m) and 
?/v(i?-m-/f), i.e., t a f tl = J d 3 #Vv(£-m)*(-^ + 

V(R))tlj cr (R — (m + If)). From the exact band struc- 
ture e^g, we found that the tunneling amplitude between 
the nearest neighbor sites along any direction is always 
much larger than other ones. For example, for Vl = 4i?R, 
Vs = H-E^, the one for the next nearest neighbor sites 
tsidxAtpidxj) are only 5%(2%) of t Syd& ^{t p4& ^) respec- 
tively. We thus keep only the tunneling between the near- 
est neighbor sites in the Hamiltonian. On the other hand, 
we have verified that introducing such a small correction 
from the tunneling between the next nearest neighbor 
sites does not change the qualitative feature of the phase 
diagram. 

There are also interaction terms for particles at dif- 
ferent sites besides on-site ones. These terms, however, 
depend on overlap integrals of Wannier wave functions 
at different sites, either in the same or different bands. 
Those integrals are always much weaker than the onsite 
interaction. For example, we have verified from exact 
numerical results that the largest nearest neighbor inter- 
action is one or two orders of magnitude smaller than the 
onsite interaction for the parameters regime we consid- 
ered. Thus we kept only onsite interaction in the Hamil- 
tonian. 
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